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Abstract 

Using the nonrelativistic effective field theory vNRQCD, we determine the contribution to the 
next-to-leading logarithmic (NLL) running of the effective quark-antiquark potential at order v 
(l/m|k|) from diagrams with one potential and two ultrasoft loops, v being the velocity of the 
quarks in the cm. frame. The results are numerically important and complete the description of 
ultrasoft next-to-next-to-leading logarithmic (NNLL) order effects in heavy quark pair production 
and annihilation close to threshold. 
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I. INTRODUCTION 



The measurement of the hne-shape of the total top-antitop quark production cross section 
in the threshold 2mt is one of the major tasks of the top quark physics program 

at a future linear collider. The most prominent quantity to be measured is the top quark 
mass, and one can expect an improvement in precision of mt by about an order of magnitude 
to mass measurements based on reconstruction obtained at the Tevatron and the LHC 
In addition, fitting the lineshape measurements to theoretical higher order predictions allows 
to control precisely the top mass scheme, which is not the case for mass reconstruction at 
hadron colliders based on Monte-Carlo generators. 

To obtain a meaningful theoretical description of the nonrelativistic threshold dynamics 
it is required to systematically sum the so-called Coulomb singular terms oc {ag/vY in a 
systematic nonrelativistic expansion, v being the relative velocity of the top quarks. This 
task is achieved by means of effective theories based on nonrelativistic QCD (NRQCD) 0]. 
In this approach, however, sizeable logarithmic terms oc (a^lnw)"' are not systematically 
accounted for, which leads to rather large normalization uncertainties of the cross section 
line-shape. The next-to-next-to-leading order (NNLO) predictions in this "fixed-order" ap- 
proach were estimated to have a normalization uncertainty of order 20% [3|. Presently 
known NNNLO fixed-order corrections, e.g. 043], seem to reduce this uncertainty to about 
10% jl]. This normalization uncertainty might not affect the top mass measurement, which 
primarily depends on the cm. energy where the cross section rises, but it still renders preci- 
sion measurements of other quantities such as the top total width or the top quark couplings 
impossible. On the other hand, such large normalization uncertainties might cast general 
doubts on the reliability of the theoretical method itself. To match the statistical uncer- 
tainties, that are expected for these quantities at an International Linear Collider (ILC), a 
theoretical precision of the cross section normalization of at least 3% would be required. 

In Refs. (see also Ref. [13]) it was demonstrated that the summation of logarithmic 
(osln^;)" terms, using renormalization-group-improved (RGI) perturbation theory, can sig- 
nificantly reduce the normalization uncertainties of the threshold cross section. Concerning 
QCD effects the RGI leading-logarithmic (LL) and NLL order predictions of the threshold 
cross section are completely known, but no full NNLL order prediction exists at present. 
The full NNLL order prediction is, however, required to obtain a reliable estimate of the 
remaining theoretical normalization uncertainties. 

The missing ingredient is the NNLL running of the Wilson coefficient of the leading order 
effective current that describes production and annihilation of a nonrelativistic ti pair in 



a S-wave spin-triplet state ('^S'l). Adopting the label notation from vNRQCD |13[| the 
current has the form 

Ji,P = V'pCr(^c^2)X-p, (1) 

where V'p and Xp annihilate top and antitop quarks with three-momentum p, respectively, 
and where color indices have been suppressed. The current does not have a LL anomalous 
dimension because there is no one-loop vertex diagram in the effective theory that contains 
UV divergences associated with the current Ji p. Such UV divergences arise at NLL order 
from insertions of the next-to-next-to-leading kinetic energy operators and from insertions 
of the next-to-next-to- leading order potentials The corresponding computations were 



carried out in Refs. 12|-|l4l] and are completed |l3l. Il4l|. Using the conventions from [I3j the 



resulting NLL order renormalization group equation (RGE) for the Wilson coefficient ci of 
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the current Ji^p has the form (S^ = 2) 



' d 



NLL 



167r2 



V. 



is). 



is) f 

k.cSi 



(2) 



where u is the vNRQCD velocity renormahzation parameter that is conveniently used to 
parametrize the correlation between soft and ultrasoft dynamical scales within the renor- 
malized effective theory The generalization of Eq. for currents describing pairs of 
quarks and colored scalars in any angular momentum and spin state (^'^^^L j) were derived 
in Ref. 

In Eq. ([2]) the term Vc^^ is the Wilson coefficient of the Coulomb potential oc 1/k^, and 
Vi'^ and V^"^ are the coefficients of the O(u^) potentials with the momentum structure 
and (p^ + p'^)/(2m^k^), respectively, m being the heavy quark mass. The term vi**'* is 
the coefficient of the spin-dependent potential that can contribute for spin triplet S-wave 
states. The coefficient V^^^^ combines Wilson coefficients from several non-Abelian l/m|k| 
potentials at 0{v). At LL order it has the form 
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The first term arises from soft renormahzation and the second represents the evolution from 
ultrasoft UV divergences.^ The terms in Eq. ([3]) are defined such that the ultrasoft terms 
have zero matching condition at z/ = 1, (Vf/|!g(l)).us = 0. The superscripts (s) refer to the 
color singlet state of the quark pair. 

At NNLL order there are two types of contributions that have to be considered for the 
evolution of Ci. The first arises from three- loop vertex diagrams that come from insertions of 
subleading soft matrix element corrections to the potentials and from insertions of potentials 
with additional exchange of ultrasoft gluons. The corresponding results were determined in 
Ref. [TgI and are referred to as the non-mixing contributions as they affect the evolution 
of Ci directly through UV-divergences. The second type of contributions arises from the 
subleading evolution of the potential Wilson coefficients that appear in the NLL order RGE 
shown in Eq. They are referred to as the mixing contributions as they affect the evolution 
of Ci indirectly. Except for the coefficient of the Coulomb potential Vc^"* [HI, 17 1 and for the 
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spin-dependent potential Vs no complete determination for the subleading evolution 

exists at present. 



The analysis of the three-loop (non-mixing) terms in Ref. [16[ showed that the contribu- 
tions involving the exchange of ultrasoft gluons are more than an order of magnitude larger 
than those arising from soft matrix element insertions and in fact similar in size to the pre- 
viously known NLL contributions. The reason is related to the larger size of the ultrasoft 
coupling as{mv'^) and to a rather large coefficient multiplying the ultrasoft contributions. 
These large ultrasoft contributions are responsible for an uncertainty in the normalization of 
the most up-to-date top pair threshold cross section prediction of at best 6% 19|, which is 
quite far from the required precision (see also Ref. [lO|). From this analysis it is reasonable 
to assume that the ultrasoft effects, which form a gauge-invariant subset, also dominate the 



^ There is a misprint in Eq.(4) of Ref. [3l where a factor 1/2 is missing for the soft contributions. All 
numerical computations in Ref. [3] were carried out with the correct result. 
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mixing contributions. This is also consistent with the small numerical effects |20| of the 
NLL evolution of the spin-dependent coefficient vi**^ [l8|, which is dominated by soft effects 
and receives ultrasoft contributions only indirectly through mixing. 



In Ref. [2l| we have computed the ultrasoft contribution to the NLL anomalous dimen- 
sions of the spin-independent 1/m? potential coefficients V^^ and vl^\ Here UV-divergent 
diagrams with two ultrasoft loops contribute directly to their anomalous dimensions. These 
NLL ultrasoft corrections were indeed large and found to be similar in size to the previously 
known LL order evolution. Their overall contribution in the NNLL mixing corrections to the 
NLL RGE in Eq. ([2]), however, turned out to be quite small due to the l/lGvr^ suppression 
factor. 

In this work we complete the determination of the NNLL ultrasoft corrections to the 
anomalous dimension of Ci by the computation of the NLL ultrasoft renormalization group 
(KG) evolution of the coefficients (V^^^fj)?!^. Since its contribution in Eq. ([2]) is substantially 

(s) 

larger than those of the V2ri the impact of the NLL ultrasoft corrections can be expected 
to be sizeable and might compensate the large ultrasoft NNLL non-mixing corrections to 



ci determined in Ref. [jJl. As we show in this work, this is indeed the CclSG clS the results 
of this work substantially stabilize the RG evolution of Ci. Our result will contribute to a 
reduction of the theoretical uncertainty of current RGI predictions for the heavy quark pair 
threshold production rate and in particular for top quark pair production at a future linear 
collider. 

This work is organized as follows: In Sec. |TT]we briefly review the effective theory setup 
used for our work and define the relevant operators in the vNRQCD Lagragian. The ultra- 
soft renormalization procedure with regard to the l/m|k| potentials is explained in Sec. IIIII 
Section |IV] gives a short overview of the respective calculation at the LL level in Feynman 
gauge, before we present the new calculation at NLL order in Sec. |Vl In Sec. I VI I we sum- 
marize the results for the l/m|k| and potentials and present the complete ultrasoft 
NNLL mixing contribution to the running of the current coefficient Ci. Sec. IVIII contains a 
numerical analysis of our results and Sec. IVIII our conclusion. 



II. THEORETICAL SETUP 



The vNRQCD Lagrangian is organized as an expansion in the heavy quark velocity v. It 
consists of three mayor parts ll|, ll3|, |22| , 



vNRQCD 



(4) 



containing kinetic terms and ultrasoft interactions potential interactions [Cp) and 

interactions involving soft degrees of freedom respectively. The ultrasoft term has 
the form 



2m 



P 



where = + iguA^{x) is the ultrasoft gauge-covariant derivative, gu is the ultrasoft 
coupling constant, and G'^'^ is the ultrasoft field strength tensor. Besides the propagation of 
the heavy quarks describes their interaction with ultrasoft gluons. 
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The potential term £p contains the potential interactions between the quark and the 
antiquark and can be written as 



pk ' 



where 



(6) 



(7) 



consists of potential four-quark operators. They depend on the soft three-momentum labels 
p and p' of the heavy quark fields, and the coefficients have the form 
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the terms a, /3, A, r being color indices and k = p' — p. The ellipses in Eq. (|8]) refer to 
spin-dependent 0{v^) potentials, that are not relevant in this work, as well as to higher 
orders in the v expansion. The leading order terms in Eq. (|8]) in the velocity power counting 

are the Coulomb potential operators. Their LL coefficients are (miy) and 

Vc^\i') = 0, where u is the vNRQCD velocity renormalization scaling parameter fll| . 

The term Cpu includes the higher order terms in the multipole expansion related to the 
potentials in Eq. ([8]). It is fixed by reparametrization and (ultrasoft) gauge invariance of 
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The first term in Eq. describes the coupling of an ultrasoft gluon to the Coulomb poten- 
tial, the other terms are four-quark operators with ultrasoft derivatives acting 
on the fermion fields to the left and to the right. The associated Feynman rules can be 
found in App. [Aj 



In contrast to early works [11|, Il2|, |22| , we omit potential terms of the form 
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in Eq. ([8]). In the term 
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we instead introduce the sum operators Oki as they are suitable to carry out the ul- 
trasoft renormahzation procedure 13|, |l6|. These sum operators generate a potential 
Vfc oc l/m(k^)^~ upon summation of the intermediate soft label momentum q, when inserted 
into vNRQCD matrix elements. The ellipses in Eq. ( ITTl) indicate other sum operators 13, 16[ 
which we, however, do not need to consider in this work. For the sum operators needed in 
our work we use the following operator basis, fis = being the soft renormahzation scale: 
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The functions fi depend on the external soft three-momentum labels p, p' and the interme- 
diate soft three-momentum label q. They are defined as 13 
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In four-quark matrix elements the sum over q is understood to be replaced by a soft loop 



integral [11 



I.e. 
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The q integrals over the /j's in d — 1 spatial dimensions yield functions that only depend on 
k = p' - p. 
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They reproduce for c? — 4 the l/|k| potentials mentioned in Eq. ( ITOj) . For simphcity we 

call the sum operators Ol^^ also Vk potentials in the rest of this work. The operators O^j}^ 

and have been introduced before in Ref. [3] for the LL ultrasoft renormalization. The 

operators O^g' are new and required for the presentation of our computations at NLL 
order. 

Unlike the corresponding operators in Eq. fllOp the sum operators can be renormalized 
consistently beyond one-loop level 13, 23|. By including the renormalized factor [vi"^^(i^)]^ 



m 



the definition of the O^i we anticipate a factorization of the soft (LL) contributions (encoded 



in Vc'^\h')) and the ultrasoft contributions (encoded in V],f ''(z/)) to their RG running. An 



analogous factorization occurs for the potentials Vr and V2 13|, |2l|.^ In this paper 



we compute the (ultrasoft) NLL running of the Vki{i^)- The LL running was determined in 



Ref. 13 



in Coulomb gauge. 



Using the equation 
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(17) 



for vectors in color space, we can easily convert the potential operators in Eqs. (E]), (E]), ( 1T2|) 
from the basis formed by the two 3 ® 3 color structures 1 ® 1 and ® into the color 
singlet/octet basis used for physical applications. In this paper only the gauge invariant 
color singlet configuration of the potentials (and the Oki) will be considered. In order to 
provide compact results we will therefore restrict ourselves to the renormalization of the 
singlet operators 
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Singlet Yl + /l + 2/2) [V^p/ Xlp' X-p] 

p,p',q 



(18) 



p,p ,q 



1^1 — 2 ^)1 ® 1 — cTx"^^ ® '^^ color singlet projection operator. 



Here Psingiet = (1 - 2 g^)l ® 1 - ® f ^ 

All fields, couplings and Wilson coefficients in the above Lagrangian are to be understood 
as bare quantities unless stated otherwise. For the renormalized quantities, indicated here 



^ A complete account on how this factorization arises requires the introduction of additional sum operators 
with four heavy quark and two soft fields, which, analogous to the Oki, absorb the UV divergences of 
vNRQCD diagrams with external soft fields and ultrasoft and potential loops. By forming a soft tadpole 
these operators contribute soft mixing terms to the anomalous dimension of the O^i- Together with the 
contributions in this work they produce the factorized running of the coefiicient [vP {'^)? M- This 



approach has been adopted in Ref. [21| . We do not follow this approach here and refer to Ref. [23| for a 
detailed discussion. In the two-stage matching approach adopted in "potential" NRQCD (pNRQCD) |24l |. 
this kind of factorization in the potential coefficients is implemented by construction. 
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by the index R, we chose the usual conventions in c? = 4 — 2e dimensions: 



9u — t^u^gdu, 9s — f^s^9 9s ■ 

= 4/^71^, Za = 1+SZa, ^^^^ 



where /i^ = mi/ and fiu = mv"^ are the soft and ultrasoft renormahzation scales, respectively. 
For convenience, we will drop the index R throughout this paper and only deal with MS 
renormalized quantities in the following. We will moreover use the notation 

9"^ ^2 
as = as{niv) = , au = asimu'^) = . (20) 
47r 47r 



III. RENORMALIZATION PROCEDURE 

The aim of this work is to determine the ultrasoft anomalous dimension of the sum oper- 
ator coefficients V^/ (u) at NLL order. We work in Feynman gauge and use the MS renor- 
mahzation scheme in combination with dimensional regularization {d = 4— 2e). All operators 

in Eq. (JH]) have zero matching condition at the hard scale, i.e. v£'"^''(i/ = 1) = (l3| . 
Purely ultrasoft loop corrections to potential operators in vNRQCD are inevitably sup- 
pressed by two powers of the velocity v.^ We therefore need at least one additional compen- 
sating factor as/v from a potential loop in order to obtain the correct w-scaling of a Vk 
potential [H, 22|. This is why the Vk potential first receives a running at as indicated 



by the factor (Vc^'')^ in the definition of the 0^^f"\ 

FigureEJand table|T]show relevant UV-divergent and diagrams, which 

contribute in Feynman gauge to the anomalous dimension of the Vk potentials at ultrasoft 
LL and NLL level, respectively. Each diagram consists of one UV-finite potential loop and 
UV-divergent ultrasoft gluon attachments. Many diagrams also contain ultrasoft derivative 
operators from higher orders in the multipole expansion. They are depicted as nabla (V) 
and Laplace (A) operator symbols on heavy quark lines or on potential vertices. The 
operator insertions on the lines represent relativistic corrections to the heavy quark kinetic 

2 

term oc (9o — ^) (see Eq. (j5])) due to ultrasoft momentum components fiowing through the 
heavy quark propagators. The potentials with ultrasoft derivatives are defined in Eq. ([9]) and 
represent corrections due to ultrasoft momentum exchange between quark and antiquark. 
The corresponding Feynman rules are given in App. In Fig. [T]we show our generic choice 
for the external energy and momenta running through the diagrams with a potential loop. 
Note that vNRQCD as reviewed in Sec. [Ill is formulated in the center of mass frame. The 
little arrows on the heavy quark lines denote positive energy fiow and are suppressed in the 
other graphs shown in this paper. 



At 0(w°) the interaction between ultrasoft gluons and heavy quarks can be removed by a field 



11 



23 



redefinition that leaves physical predictions unchanged 
^ Such derivative operators also exist for heavy quark lines /vert ices which do not carry ultrasoft momentum 
components, but then they evaluate to zero. This needs to be taken into account when showing the 
equivalence of diagrams that differ by the routing of ultrasoft momenta. 
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FIG. 1: External energy/momenta convention for all diagrams with potential loop. 



As demonstrated in Ref. 13[ the additional potential loop entails the renormalization of 



the Olj^ sum operators instead of the potentials in Eq. ( ITOl) . The reason becomes apparent 
at the three-loop level, when the potential loop appears together with two ultrasoft loops. 
Some typical example diagrams are shown in table HI After the integrations over ultrasoft 
and potential loop momenta their leading divergences typically take the schematic form [l^ 

Ikl 



k\\k^) 2e^[E^) e^\E^) ~ \k\[ 2e^ 2e ^""Ik^ + ' ' ' J ' ^^^^ 



where E represents the ultrasoft scale ~ mv'^, i.e. the kinetic energy of the external heavy 
quarks. Inside the square brackets the first term arises from the three-loop diagrams and 
the second from two-loop diagrams with one-loop counterterms to remove sub divergences. 
The factor in front of the square brackets on the left-hand-side comes from the finite po- 
tential loop. The ln(/i|/k^)-dependent 1/e-pole on the right- hand- side of Eq. fl2T]) cannot 
be absorbed into counterterms for potentials of the type in Eq. ffTOj) . because this would 
cause inconsistent anomalous dimensions. This is the reason, why the sum operators Oki 
are essential for a consistent renormalization beyond the one-loop level in vNRQCD. Upon 
summing the intermediate soft momentum label q, they contain the proper potential pref- 
actors in Eq. (^T^ in their operator structure. Therefore the 6Vki counterterms only absorb 
proper, scale independent UV poles and lead to consistent anomalous dimensions. 

In practice we have to determine the 6Vki from the ultrasoft UV divergences of two- (LL) 
and three-loop (NLL) Feynman diagrams after the integration over the zero-component of 
the potential loop momentum (go), which can be carried out using residues. The potential 
three-momentum (q) integration, however, is not carried out in the renormalization proce- 
dure as the q-dependence of the remaining UV-divergences has to be matched onto the label 
structures of the sum operators Oki- There is a technical subtlety concerning the vNRQCD 
multipole expansion related to taking residues with respect to the potential loop momen- 
tum go in the upper or the lower complex half plane. Here differences can arise in the soft 
momentum label structure in terms of p, p' and q resulting from insertions of the ultrasoft 
derivative operators which are non-zero only when the corresponding quark line carries an 
ultrasoft momentum component. Upon carrying out the label sum the results are unique. 
To avoid redundancies in the operator structure of the Oki we need to define a consistent 
prescription for the respective Feynman graphs in order to ensure, that the operator basis 
of the Oki required for the renormalization is unique. We evaluate each Feynman diagram 
with a potential loop considered in this paper as follows (g'^ = (go, q) being the potential 
loop momentum): 

1. Neglect the kinetic and potential multipole correction operator insertions for the mo- 
ment. Pick one internal heavy quark or antiquark line (which carries momentum q) 
and choose an ultrasoft momentum routing such that this heavy quark/antiquark line 
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does not carry any ultrasoft momentum. Now, according to the chosen routing insert 
the required ultrasoft derivative operators in all possible ways. 

2. Compute the residue of the pole in the complex go plane associated with the propa- 
gator without ultrasoft momentum picked in step 1 for the chosen routing and each 
vertex/insertion configuration. Multiply with the proper factor {±27ri) according to 
the residue theorem depending on whether the pole lies in the upper or lower complex 
half plane. 

3. Perform the ultrasoft loop integrals. 

4. Repeat steps 1-3 for all internal heavy quark/ ant iquark lines. Sum all expressions 
obtained in this way and divide the result by two. 

The factor 1/2 from point 4 compensates for the overcounting, which arises in points 1 and 2 
from closing the go integration contour in the upper as well as in the lower complex half 
plane. 

From a technical point of view this procedure corresponds to the renormalization of 
properly defined six-quark operators in order to absorb the ultrasoft UV divergences of the 
six-leg diagrams, which are generated from cutting the heavy quark/antiquark lines of the 
potential loop where the residues of step 2 are taken. This is the origin of the notation used 
in Fig. |2l The two heavy quark fields associated with the cut can then be contracted to form 
tadpole diagrams with four external legs, which is equivalent to carrying out the label sum 
for the sum operators. This procedure yields the same results as using the matrix elements 
of the sum-operators (see Ref. [1^ for details). In Fig. [2] we visualize the terms arising in 
point 2 of the above caculational prescription using their analogy to the described six-leg 
diagrams. 

IV. LL CALCULATION 




FIG. 2: Example for a calculation with one ultrasoft loop involving ultrasoft derivative operators. 
Each of the cut heavy quark/antiquark lines on the right-hand side corresponds to a qo residue to be 
taken at this line within the original closed potential loop diagram on the left-hand side according 
to step 2 of our calculational prescription. The color structure of the diagrams is understood to be 
not affected by the cutting, i.e. it is the same on both sides of the equation. 
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In order to illustrate our computational procedure let us consider the simple case of the 
two 0{a^au) diagrams shown in Fig. [2l In Feynman gauge they contribute to the LL 
anomalous dimension of the Oki and contain ultrasoft derivative operators. The outcome 
of steps 1 to 4 are sketched graphically in Fig. |2J Each six-leg diagram corresponds to 
picking up the residue of the cut heavy quark/antiquark line as described in step 2. The 
overall factor 1/2 is due to step 4. Including in addition all possible mirror-graphs and the 
analogous diagrams with one A-operator, we obtain (up to an overall color factor not shown) 

the result + C'(e°), where in the integrand the term (/o + /i)/4 arises from 

the residue in the upper complex half plane and (— 3/o + 5/i)/12 from the residues in the 
the lower complex half plane. Upon integration over q both residue contributions agree. 

We now briefly discuss the full LL ultrasoft renormalization of the Oki in Feynman gauge; 
for details see Ref. jlH. This calculation is also interesting by itself since the available 
vNRQCD computations of Refs. [13|, |22| were carried out in Coulomb gauge. Fig. [3] shows 
the relevant diagram topologies. Possible insertions of ultrasoft derivative operators are not 
shown. 




FIG. 3: Diagram topologies that generate the RG running of the Oki at LL level. Up-down and 
left-right mirror graphs are understood. 



The ultrasoft gluons in diagrams E] a-c are spatial (A type) because there are no four-quark 
operators with one additional ultrasoft gluon in vNRQCD. Thus the vertices in these 
graphs already contribute two powers of v to the velocity counting and we do not need to 
insert further derivative operators. On the other hand, diagrams with a spatial A gluon 
and topologies d-j are either UV finite or their divergence is exactly canceled by one loop 
diagrams with 6Vr, SV2 or wavefunction counterterm insertions. The latter are depicted in 
diagrams [3] k and 1, respectively. The very same applies to the diagrams with an gluon for 
topologies d-j, where two kinetic V-operators are inserted on the heavy quark lines within 
the ultrasoft loop. 

The remaining diagrams d-j with an A^ gluon having the correct v scaling contain either 
one A- or two V-operator insertions, of which at least one is the V-potential in Eq. (lA5p . 
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Adding them to the A diagrams |3] a-c gives 
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(22) 



The result is exactly the same as obtained in Refs. |l3l . |22[ where Coulomb gauge has been 
used. We note that the Coulomb gauge computation at LL order is substantially simpler 
since in Coulomb gauge gluons cannot propagate and the number of diagrams that have 
to be considered is reduced considerably. However, the Coulomb gauge calculation becomes 
a nightmare at the two-loop level and seems entirely unfeasible. We also note that we 
regularize IR divergences by adopting an off-shell configuration for the external heavy quark 
lines with E ^ p^/2m = (p')^/2m. 

From Eq. (!22|) we can easily read off the already known 0{au) color singlet counterterms 
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V. NLL CALCULATION 

Let us now consider the divergent 0{a'^g a^) diagrams contributing to the ultrasoft NLL 
anomalous dimension of the V^j, i.e. diagrams with one potential and two ultrasoft loops. 
In table [T] we have classified the relevant diagrams without gluon selfenergy with respect to 
their vertices and kinematic operator insertions. For each class one sample diagram with 
a specific configuration of vertices and insertions is shown graphically. Diagrams with all 
possible permutations of the vertices and different attachments of the gluons to the heavy 
quark lines are understood to belong to each class. The respective contribution in the last 
column of table [T] denotes the sum of the divergent parts of all class members. Contribution 
I covers classes 1 and 2, contribution iv covers classes 12 and 13. In the following we will 
call the diagram classes 1-8 "Abelian" and 9-14 "non-Abelian" . The relevant Feynman 
gauge Feynman rules are found in Appendix |Al 

Besides the Abelian and non-Abelian diagrams in table [T] there are (!?(a|Q;^) diagrams 
containing one-loop gluon selfenergy subdiagrams with ultrasoft light fermion, ghost and 
gluon bubbles. For the calculation of these diagrams we take the graphs with only one 
ultrasoft loop in Fig. [3] and replace the ultrasoft gluon propagators by their one-loop cor- 
rection. Note that in Feynman gauge the selfenergy insertion is a non-diagonal matrix in 
the components of the ultrasoft gluon field and can couple the zero-component {A^) to the 
spatial components (A) of the ultrasoft gluon field due to its relativistic nature. Therefore 
we have to take into account diagrams with one A^-vertex, one A-vertex and a single V- 
operator insertion, in addition to the diagrams we obtain by simply inserting the selfenergy 
bubble into the 0{a\au) diagrams of the previous section. After subtraction of the selfen- 
ergy sub divergences by adequate counterdiagrams the sum of all ultrasoft UV divergences 
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from the three-loop a^) diagrams with a gluon selfenergy, reads: 



Ca{5Ca - 4:Tnf) Ca{31Ca - 20Tn 



f) 



72e2 216e 



m TT^ 



Ci.(CA-2C^)(/o + /i + 2/2)l®l + (C^-4Ci.)(/i + /2)T^®T^ , (24) 



(27r)'^-i 

where T = | in the fundamental representation of S'f/(3) and n/ is the number of light 

quark flavors. The result contributes to the counterterms 5Vl^ and 5V^2 • The gauge 
invariant contributions from light fermion loops oc n/ in Eq. (!24|) were also determined in 



Coulomb gauge [25 



In table [TTl we present results for the UV divergent parts of the Abelian and non-Abelian 
diagrams classified in table |T]after all one- and two-loop sub divergences that contribute to the 



renormalization of the 1/m potentials [12|,|21[ have been subtracted. The net contributions 



from the Abelian diagram classes 3 and 4 vanish, because the diagrams are either finite or 
the UV-divergences are cancelled exactly by the respective diagrams with potential 
counterterms. The results for I and II-I-III+IV+V also contain contributions from external 
wavefunction diagrams, i.e. heavy quark selfenergies on the four external legs. Two-loop 
heavy quark selfenergies with a triple gluon vertex vanish. 

In order to avoid singularities of infrared origin we again assigned a uniform off-shellness 
to each of the four external legs. After the subtraction of sub divergences the logarithmic 
dependence on the off-shellness vanishes. 

Note that apart from the external off-shellness the ultrasoft loops in the six-leg diagrams 
also involve the physical off-shellness of the intermediate (uncut) heavy quark line, 2{E — 

2 

^). So the nontrivial ultrasoft two- loop diagrams that have to be calculated represent 
complicated multi-scale integrals. 

We solved the integrals using partial fraction decomposition and integration by parts 
techniques. To handle the larger number of diagrams (corresponding to (9(10^) six-leg 
diagrams) we developed Mathematica codes for the generation and the evaluation of the 
corresponding amplitudes. For the e = expansion of the hypergeometric functions 



resulting from the two-loop integrations we used the Mathematica package HypExp 1.1 [26 
and checked the results numerically. In order to automatize the determination of the color 
structures of the different three-loop topologies we moreover employed a few routines from 
the FeynCalc 4.1 package 

From the 0{afj) UV divergences in table [III and Eq. we can now determine the 



5 



The results in table |TI] can also be found in Eqs. (6.33) and (6.34) (Abelian contributions without wave- 



function renormalization contributions) and table 6.9 (non-Abelian contributions) of Ref. 23[. In table 6.8 
of Ref. [2^ a result for the non- Abelian contributions is shown that is based on a calculation of six-leg 
diagrams, as explained at the end of Sec. lIIH using an ofF-shcll configuration for the two lines coming from 
cutting the intermediate heavy quark propagator. This ofF-shcll configuration disagrees with the on-shell 
condition arising from doing the potential energy integration by the residue theorem, and the result is not 
relevant for the computation intended in this work. 
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class 


example 


diagram 


vertices 


contrib. 


potential 


quark-gluon 


3 gluon 


kin. ins. 


1 




1 X Vc 

1 X A-k 


1 X A 

2 X ^0 






I 


2 




2 X A-k 


2 X ^0 






3 






2 X 


2 X A 
2 X 









4 






2 X Vc 


4 X ^0 




2 X V-p 





5 




2 X 


4 X ^0 




1 X V2 


II 


6 




2 X V-k 


4 X ^0 






III 


7 




1 X V, 
1 X V-k 


4 X ^0 




1 X v-p 


IV 


8 




1 X 


4 X ^0 






V 


9 






2 X Vc 


2 X A 

1 X yl^ 


A2^0 




i 


10 






2xVc 


1 X A 

, , .ft 


A(A0)2 


1 X v-p 


ii 


11 




1 X Vc 

1 X V-k 


1 X A 


A(A0)2 




iii 


12 






1 X Vc 

1 X A-k 


1 X A 

1 X A^' 


A2^0 




iv 


13 




2 X A-k 


1 X A^ 


A2^0 
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1 X Vc 

1 X A-k 


2 X ^0 


A(A0)2 


1 X v-p 


V 
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1 X V-k 
1 X A-k 


2 X ^0 


A(A0)2 




vi 



TABLE I: Diagram classes with one potential and two ultrasoft loops (without gluon self energy). 
Diagrams with all possible permutations of the vertices and different attachments of the gluons to 
the heavy quark lines are understood to belong to each class. The vNRQCD Feynman rules for 
the respective vertices, propagators and kinematic operator insertions are given in App. [Al 

14 



contribution 



result -(^ 



- J^)^[CF(CA-2CF)(/0+/l+2/2)155i + (CA-4CF)(/l+/2)T'4®f'*] 



II+III+IV+V 



(~ ih^ + ) ^ [Cf (Ca-2Cf)(/o+/i+2/2)1®1 + (Ca-4Cf)(/i+/2) T^<S>f^] 



1+n 



48.(q-p)^(p^-q)^(2mj;-q^) KCC^I-SCaC^+SC^ + -6C^Ca + 12C2 )r^®f a] 



16 77 



111 



^Ca [{CA-2Cp)CF(fo+fi)l&+{CA~4:CF)hT'^m^] + 
+ 5^[Cj.(CA-2CF)l(g)l{/o[6CA+{2C'A-12C'F)7r2]-2/i[-3CA+(2C'A-6C'F)7r2]}- 
+ T-^^f { 12Cf (2CF-CA)7r2/o+/i [6Ca (CA-4CF)-4(Ci-4CFCA+6C|.)7r2] }] 



IV+V+Vl 



if^ ^C-A [(CA-2Cj-)CF(/o+/i)lCS5l+(CA-4Ci.)/ir^®f^] + 
+ 5^ [cF(CA-2CF)l®l{-5/o[6CA + (2CA-12CF)7r2]+10/i[-3CA+(2CA-6CF)vr2]-4CA/2[12-27r2]} 
+r'*(g)f'*{60{CA-2CF)CFvr2/o+/l[-30CA(CA-4CF)+20(C2-4CFCA+6C|.)7r2] + 
+8CA/2[-3CA+12CF+{2CA-2CF)7r2]}] 



TABLE II: UV divergent contributions before the q integration from three-loop diagrams with one 
potential and two ultrasoft loops as defined in table [H Subdivergences have been subtracted by 
one- and two- loop counter diagrams. The soft three-momentum structure is (mainly) encoded in 
the fi functions of Eq. ([13]) . 



sum-operator counterterms SVki- They take the generic form 



(25) 



In addition to Eq. ( !23l) we included here, following Refs. [13| and [28[, so-called soft pull-up 
terms at 0{as) and 0{ag). They originate from the zerobin subtractions in corresponding 
amplitudes with one and two soft loops j28|. The ultrasoft one- loop coefficients Ai were 
given in Ref. [lH and can be read off Eq. ( 12^ . They read in the color singlet channel 



C aC f{C A — '^C f) 



A 



CaCf{Ca—^Cf) 



37r ' " Svr ' 

For the anomalous dimension of the Vki{v) to be finite, the relation 

/3o 



A 



is) 



0. 



Svr 



Ai 



(26) 



(27) 



is a necessary condition. 
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FIG. 4: Mercedes star type diagram contained in diagram class 9 (in table |T]). Together with the 
Mercedes star diagrams of class 10 it gives rise to the non-Abelian contribution (i+ii) in table [III 
The upside down graph is understood. 



Summing all Abelian and non-Abelian contributions in table HI] except for contribution 
(i+ii), which we will discuss later, and projecting the result on the physical color singlet 
channel we find: 



I + II + III 



X 



IV + V + iii + iv 



VI 



^2 



1 



7 

2Ae 



7T 

18e 



(T)^2..4. 



Cp(Ca-2C^)(/o+/i+2/2; 



12e2 24e i»e rmr^ 
-CFiCA-^CF)if, + h) 



Singlet 



(28) 



The C(a^/e^) terms of Eq. ( l24l) and Eq. ( l28l) satisfy Eq. ( 1271) . which provides an impor- 
tant cross-check of our calculation. From the respective (9(a^/e) terms we determine the 
coefficients 



B. 



is) ^ Ca{Ca - 2Cf)Cf[Ca{^7 + Qn^) - lOnjT] 
^ 1087r2 
(.) _ CAjCA - ACf)Cf[Ca{^7 + 67r2) - 10 njT] 
^ 1087r2 



(29) 



The contribution (i+ii) in table [TTl is generated exclusively by those diagrams of class 9 and 
10 which belong to the subclass with the Mercedes star topology shown in Fig. HJ All other 
diagrams of class 9 and 10 add up to zero. The result requires a separate discussion. It has 
already been computed in Ref. |23|, but not considered there for numerical examinations. 
Although its effects are numerically tiny, they are interesting conceptually and worth to 
explain in some detail. Up to a color factor the term (i+ii) has the form 



48e (q 



p)2(p' 



q)2(2mE - q2 



(30) 



If one carries out the renormalization procedure based on 5-loop current-current correlators 
instead of quark- ant iquark scattering amplitudes, as realized in Ref. 16| for the computation 
of the 3- loop NNLL non- mixing contributions to the anomalous dimension of Ci, this term 
generates additional ultrasoft NLL terms in the counterterms of the l/m2 potential coeffi- 
cient Vr and of the coefficient of the l/m2-suppressed S-wave current C2. These contributions 
have not been considered in the literature before. They were also not contained in the results 
of Ref. 21[ since they do not arise for scattering diagrams with equal external off-shellness, 
which was used for regulating IR divergences in Ref. [Hj . Figure [5] shows a sketch of one 
of the 5-loop current-current correlator diagrams involved in the renormalization procedure 
following Ref. [l6| . The divergences of these 5-loop diagrams have to be canceled by 3-loop 
current- current diagams without ultrasoft gluons, but with one of the potentials/currents 
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FIG. 5: 5-loop current-current correlator with a Mercedes star subdiagram. It is relevant in 
the current-current correlator based renormalization procedure associated with the divergent term 
(i-|-ii) in table [III In this approach one also has to take into account the 5-loop diagrams, where 
the Mercedes star is inserted into the left and right potential loop. The latter (partly) contribute 
to the renormalization of the 1/m^-suppressed production/annihilation currents. 



replaced by the respective 1 / suppressed counterterm. This approach yields an additional 
contribution to the counterterm 5Vr*'' to be added to the expression given in Ref. 21|. It 



results from insertions corresponding to the term (i+ii) in table [IT] into a current-current 
correlator and reads 



Alternatively one might write Eq. ( 130|) as 



ASVl'^ = 2(4vr«5) ASW ; Ai?^ = . (31) 



24 

/o + /i , 2mE 



96 e ' 48e(q-p)2(p'-q)2(2mE-q2) '^^^'^ 

and directly absorb the ultrasoft divergence oc /o + /i into the counterterms 6Vk3, while the 
counterterm SVr^^ of Ref. 21 1 and all other results in that reference remain unchanged. This 
leads to the coefficient 

5^ = ^. (33) 

The second term of Eq. (l32i) cannot contribute to the anomalous dimension of Ci , because the 
integrations over the soft momenta q, p, p' give finite results, when the structure is inserted 
into a current-current correlator. It would rather contribute to the running of C2, which 
becomes clear by rewriting its numerator as 2mE = l/2(p^ + p'^) plus scaleless integrals 
(over p and p'). The remaining divergence could then be absorbed by the counterterms Sc2 
of the p^/m^-currents to the left and to the right of the current correlator. 

The ambiguity in the renormalization procedure for the contribution (i+ii) is related 
to the freedom of performing field redefinitions and does not affect physical observables 
like the running of Ci. The contribution from (i+ii) to the Ci evolution equation Eq. ([2]) 
is independent of whether we use the counterterm 6Vr or 6Vk3 to subtract the relevant 
divergences. As we will see below. Bp and -B3 enter the final expression for the NNLL 
anomalous dimension of Ci in the combination CpBp + (see footnote [7]). Thus both 
approaches lead to the same result. Together with Eq. ( 12^ . Eq. fl^ and fl5^ represent the 
main results of this work. 



VI. SUMMARY OF RESULTS 



From Eqs. (ITU]) . (12 5 p and fl27p we obtain the anomalous dimensions 

""^^'^ " " "li^^^''' " ^ ^^"^ -as) + 4: Bi {2al - «|) . (34) 
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To integrate Eq. we can assume NLL running of au and as in the Ai and LL running 
in the Bi terms. Respecting the matching conditions Vki{v= 1) = 0, the solution of Eq. (|3^ 
at NLL order precision reads 



■>;NLL^ 



~A.^r.^+(%A-^-lB^[au-as), (35) 



where /3o = ^-CA — ^Tn 



a _ 34^2 
Pi - -3^A 



ACpTuf — ^CaTui are the standard one- and two- 



loop coefficients of the QCD beta function. Eq. (1351) is simple, because it only accounts for 
purely ultrasoft effects (plus the corresponding pull-up terms). We recall that the coefficients 
V^^^ contribute to the NLL anomalous dimension of Ci with an additional factor [Vi"^^'^^(z/)]^ 



according to the definition of the Oki sum operators. 

The corresponding expressions for the 1 / potentials are given by j2l| 



(vr^(.)) 



NLL 

us 

NLL 

us 



LL/ 



LL/ 



(36) 
(37) 



In this work the soft Coulomb factor vi^^'^^(z/) has not been included in the poten- 
tial operator structure in Eq. ([H]), but appears directly in the Wilson coefficients V2,r in 
Eqs. (136 p and fl37p . The square brackets have the same form as Eq. (13 5 p and represent the 
ultrasoft contributions (plus the corresponding pull-up terms) up to NLL order. In addition 
to Eqs. fl36|l and fl37|) there are also purely soft contributions to V2,r- They are fully known 
at LL order and can be found in Ref. (l2| . 

Adopting the convention in Eq. (l3T]l . i.e. adding the extra contribution ABp^ to 



the singlet coefficient Bp^ determined in Ref. [21|'' instead of introducing an additional 



coefficient we can summarize the coefficients in Eqs. ( 135|) -( 137|) as follows: 



where 



As 

b[' 

A' 

Bi' 



AM 
B^ 



A^' 
B^ 



C pic A — 2C^ 
—CaCf 
—CaCf{Ca — 2Ci? 
CaCf{Ca — 4:Cf 



AM 

bm 

AM 

bm 



(38) 



37r IOStt^ 



(39) 



^ For the results in Ref. 



21 



the scheme choice of Eq. ([55)1 was employed. 
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This simple form is quite suggestive because the ultrasoft corrections to the color singlet 
Wilson coefficients of the l/m|k| and the potentials are universal up to color factors. 

This appears to be a quite non-trivial fact given the complexity and differences in the 
computations for the l/m? and l/m|k| potentials. We emphasize again, however, that the 
form of Eq. fl38p is depending on the scheme choice associated with Eq. (131 p . 

Using Eqs. (ITU]) and f llSp we can now write down the final expression for the ultrasoft 
contributions (including soft mixing and pull-up terms) to the effective color singlet Vk 
potential in Eq. ([3]) up to NLL order. The Wilson coefficient 



(vSff(-)) 



NLL 



87r2 



V) 



—CaCp{Ca + 2C^ 



[V, 



kl 

(T),LL/ xn2 



V 



NLL 



(40) 



87r2 



-^AWln-^ + 

Po OiS 



PI 



in 



/So 



appears in the NLL anomalous dimension of the Wilson coefficient Ci(z/) of the ^5"! heavy 
quark current in Eq. (|2]). The LL running of the Coulomb coefficient Vj(i^) is purely soft 
and given by Vc''"^'^^{u) = Anas- Together with the ultrasoft NLL running of the 1/m^ 
potential coefficients V2,r in Eqs. f l36|) and fl37|) . Eq. fHOj) constitutes the NNLL ultrasoft 
mixing contributions to the anomalous dimension of Ci. 



It is now straightforward to derive from Eq. (|2]) the complete two-loop ultrasoft part of 
the NNLL mixing contributions to the running of the ^5*1 current coefficient Ci. We can 



parametrize the form of Ci as [16 



In 



Lci(l) 



(41) 



where ^^^^ is the well known NLL order contribution 13|, |lj] and ■Cnm^^ NNLL 
non-mixing contribution determined in Ref. 16|. The ultrasoft part of the NNLL mixing 
contribution ^^^^ is new and reads 



fNNLL 
'?m,usoft 



7 

4^6 



^ 1 - In 



In z 



2Li, 



where^ 



Po 



3-2z 



41n(2 - z) 



A = Cf{A^^^ + 2^^)) + 3A['^ + 2A^^^ = -Cf{Ca + Cf){Ca + 2C^)A(^ 



B = Cf{bI^^ + 24^)) + 3B['^ + 2B^'^ = -Cf{Ca + Cf){Ca + 2Cf)B^'^ 



(42) 

(43) 
(44) 



^ If one uses the counterterm SVk3 to absorb the divergent contribution i+ii, the RHS of Eq. (|40|) reads 
-CpiCA + Cf){Ca + 2Cf)B^'1 
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with 



.^('^4^V'= (l + ^M^lnX\ (45) 



VII. NUMERICAL DISCUSSION 




FIG. 6: Renormalization evolution in of V2 + Vr and ^g. In panels a and b the ultrasoft LL 
(dashed lines) and LL+NLL evolution (solid lines) are shown assuming zero matching conditions 
at V = 1. The curves include, however, the soft LL running of the factors Vc and [Vc ]^ and 
pull-up terms as in Eqs. (f36]) . (i37l) and (fiOl) . Panels c and d show the corresponding plots with the 
pure soft LL contributions and the corresponding matching conditions added to all curves. 

In order to demonstrate the rather large size of the NLL order ultrasoft corrections to 
the anomalous dimensions we plot in Fig. [6] the evolution of the coefficients of the 
potentials (V2*'' + vi^^) and of the l/m|k| potential V^*]g at LL (dashed lines) and at NLL 
order (solid lines). These two sets of potentials are directly affected by ultrasoft effects. Since 
they appear explicitly in the NLL anomalous dimension of the ^Si current, they constitute 
the most important ultrasoft mixing effects in the anomalous dimension of Ci, see Eq. ([2]). 
In panels a and b only the LL and NLL ultrasoft running (including soft mixing and pull- 
up terms) according to Eqs. fl36l) . fl371) and fHOl) is displayed. For all curves zero matching 
conditions at z/ = 1 has been adopted. For both V^'^ + V^'^ and V^'lg the ultrasoft NLL 
contributions are rather large and can even exceed the LL ones. The situation remains in 
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Figs. [6t and d where we also include the full LL soft contributions and matching conditions 



from Ref. [l2j. Due to the small prefactor ~ Vc*''/ (IGvr 



the contribution from V9 + Vr ' is 



relatively suppressed compared to the effect of Vfc^efr in the evolution of Ci, Eq. (|2]). Therefore 
the ultrasoft NLL corrections to Vfc,efr obtained in this work have a substantially larger impact 



(s) 

on ci than those for Vo ~ obtained in Ref. 21 



The interplay of the ultrasoft NLL evolution of the potential coefficients V2, Vr and Vk,eS 
in Ci (see Eq. ([2])) and the NNLL non- mixing contributions in the running of Ci is examined 
for the case of top and bottom quark pairs in table llllt which represents an update to 



table III of Ref. 16|. The values of lnci(z/)/lnci(l) at the NLL and NNLL level are given 





m = 175 GeV 


m = 4.8 GeV 


1/ 






Snm V / ' Sm,us \ / 






tNNLL/.AitNNLL/.A 
?nm (^J + ?m,us [^) 


1.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.9 


0.0031 


0.0018 


0.0018 


0.0127 


0.0158 


0.0153 


0.8 


0.0066 


0.0040 


0.0039 


0.0267 


0.0360 


0.0335 


0.7 


0.0105 


0.0066 


0.0063 


0.0419 


0.0629 


0.0555 


0.6 


0.0150 


0.0099 


0.0091 


0.0582 


0.1008 


0.0824 


0.5 


0.0200 


0.0142 


0.0126 


0.0743 


0.1585 


0.1154 


0.4 


0.0259 


0.0202 


0.0171 


0.0852 


0.2584 


0.1530 


0.3 


0.0328 


0.0294 


0.0230 


0.0689 


0.4789 


0.1671 


0.2 


0.0401 


0.0460 


0.0313 








0.1 


0.0405 


0.0916 


0.0411 









TABLE III: Numerical values for the contributions to the running of lnci(i')/lnci(l) as defined 
in Eq. (|^T]) : '?nm^^(^) ^'^d 'Cm^s^(z^)- The values for m are pole masses. The numbers 

are obtained by evaluation of the analytic results using one- loop running for a^^ (in Eq. (j45p ) and 
taking ai"-^"^^ (175 GeV) = 0.107 and ai^-^^^^.S GeV) = 0.216 as input. 

as a function of i' for the top quark (m = 175 GeV) and the bottom quark [m = 4.8 GeV) 
cases. The numbers in the third (sixth) column show the sum of the complete (soft and 
ultrasoft) NNLL non-mixing and ultrasoft mixing contributions. The latter terms were 
computed in this work. As a comparison we display in the second (fifth) column the non- 
mixing contributions alone. The first (fourth) column gives the number at the NLL order. 
As anticipated, we see that there is a substantial cancellation between the NNLL mixing 
and non-mixing contributions, particularly for small values of z/ ~ a^, where the NRQCD 
matrix elements are being evaluated. This behavior is insofar remarkable as the NLL non- 
mixing contributions of Ci arise from 3-loop vertex diagrams while the mixing corrections are 
diagrammatically related to 4-loop vertex diagrams. This demonstrates in a quite coercive 
way the importance of properly summing logarithms and of accounting for mixing effects 
in NRQCD predictions. The relation of mixing and non-mixing contributions to 3-loop 
and 4-loop vertex diagrams, respectively, also explains why the mixing contributions are so 
much smaller than the non-mixing ones close to the matching scale i/ ~ 1 where velocity 
logarithms are small. 

A visual display of the results is given in Fig. [TJ where the evolution of ci(z/)/ci(l) is 
shown for top quark pair production close to threshold (m = 175 GeV). The blue dotted 
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line represents the full NLL curve known before from Refs. [13|, [IJ]. The red dashed line 



includes in addition the NNLL non-mixing contributions computed in Ref. [16| and exhibits 
a clear instability for i/ ~ ~ 0.15, where the low-energy matrix elements relevant for top 
threshold production have to be evaluated. This instability is the cause of the relati vely 



large 6-10% uncertainty in the normalization of the top pair threshold cross section [lO|, [19 
Finally, the black solid curve shows Ci(i^)/ci(l) including the full NLL and NNLL non-mixing 
contributions as well as the ultrasoft NNLL mixing corrections finalized in this work. The 
stabilization of the renormalization group evolution for ~ ~ 0.15 is obvious and will 
contribute to a reduction of the current uncertainties in RGI predictions for heavy quark 
pair threshold production. A detailed analysis of top and bottom quark production will be 
the subject of future publications. 



1.14 



1.12 - 
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1.08 - 



ci[v] 
c,ll] 




1.06 



1.02 - 



1.00 



FIG. 7: Renormalization group evolution of the ^Si current Wilson coefficient ci(z/) normalized 
to ci(l) for the top quark case (m = 175 GeV). The dotted (blue) line represents the full NLL 
result exp[,^'^'"'"], the dashed (red) line includes in addition the NNLL non-mixing contributions, 
expl^^'^'^ + C^m^^]- The solid (black) line accounts for the full NLL and NNLL non-mixing con- 
tribution as well as for the ultrasoft NNLL mixing corrections, exp[^'~^^^ -|- CSn^^ + ■^S^soft]- "^^^ 
typical value for the velocity scaling parameter relevant for the evaluation of top threshold matrix 
elements is ~ ~ 0.15. For the plot we have used as (175 GeV) = 0.107. 

Nevertheless, one should not leave unmentioned that the NNLL corrections are obviously 
still quite large. In particular, scale variations obtained from the NLL order result around 
ly = 0.15 fail to give a reliable estimate of the perturbative uncertainty. For heavy quark 
threshold production such behavior is not an unfamiliar situation. For example, it is well 
known that the insertion of 0{v'^) suppressed potentials and operators in NNLL order matrix 
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elements leads to large corrections that cannot be anticipated from scale variation of the 
NLL results. In the case of the renormalization scale running of ci(i^) the rather large NNLL 
shift arises from the lower order logarithms generated by the non-mixing contributions. 



VIII. CONCLUSION 

We have determined the two-loop ultrasoft corrections to the NLL order renormalization 
group equation of the 0{v) vNRQCD l/m|k| potentials. Combined with the previously 



known results for the Olv"^) potentials Vr'^ and V2\i^) determined in Ref. j2]J our result 
completes the ultrasoft NNLL mixing part of the renormalization group evolution of the 
leading ^5*1 current that describes heavy quark- ant iquark pair production at threshold. 
Together with the previously known NNLL non-mixing corrections our results present the 
parametrically dominant contributions of the NNLL anomalous dimension of the leading 
^5*1 current and should allow for meaningful RGI predictions for heavy quark threshold 
production at the NNLL level. 



Note added: While this paper was finalized we obtained a letter by Pineda [29| on the 
NLL ultrasoft renormalization group running of the spin-independent 0{l/m) and (9(l/m^) 
potentials within the pNRQCD effective field theory framework. Including the soft pull-up 
terms, see e.g. Eq. fl25p . we find agreement with our results, i.e. the previously known 
ultrasoft NLL results for the vNRQCD potentials ¥2^^ and Vr^'^ j2l| together with the 
l/m|k| potentials obtained in this work. In particular our results in Eqs. fl5^ . (^7}i and fl4Up 
agree with the corresponding position space expressions 6V}'^^Q{r; Ug, Up), 5V 2^^^g(r; i/^, t'p) 

and SV}^^Q{r; Us, Vp) of Ref. [2^, respectively, upon transformation to momentum space and 
after imposing the correlation Vp = / m = fiu{= mu'^) between the pNRQCD soft matching 
scale Ug and the pNRQCD ultrasoft renormalization scale Up. We remind the reader that 
the vNRQCD subtraction velocity u is dimensionless whereas the pNRQCD scales Up and 
Us as well as the ultrasoft vNRQCD scale fiu have the dimension of a mass. 



were more 



The two-loop ultrasoft computations in our present work and in Refs. [21|, 1231 
involved than the corresponding ones in Ref. [29| using pNRQCD. This is related to the 
simpler power counting of pNRQCD where soft fluctuations are decoupled from potential 
and ultrasoft ones. In the vNRQCD formulation used in our work soft, potential and ultrasoft 
fluctuations are contained in the theory which allows to systematically treat higher order 
mixing effects involving all three types of fluctuations (see Ref. M for such a computation). 
The equivalence of the results obtained in Ref. 29| and the ones in this work (and in 
Refs. jlll, [2^) is technically non-trivial and indicates that the structural differences in the 
vNRQCD and pNRQCD formulations might not affect the description of purely ultrasoft 
effects. 
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Appendix A: Feynman Rules of vNRQCD 



In this Appendix we list the vNRQCD momentum space Feynman rules in Feynman 
gauge, which are relevant for our calculation. Soft (±p, ±p') and ultrasoft momenta {k, ki) 
of the heavy (anti) quarks are understood to point into the direction of the positive energy 
flow indicated by little arrows on the heavy quark lines. If not explicitly indicated in the pic- 
tures the default soft and ultrasoft momentum of a heavy (anti)quark is p and /c, respectively. 



Propagators: 

ultrasoft gluon: 

ultrasoft A gluon: 
heavy (anti) quark: 

Heavy quark kinetic insertions: 



A 



k"^ + ie 

k"^ + ie 
i 



PI 

2m 



le 



(Al) 



Heavy quark - ultrasoft gluon - vertices^ 



p ■ k 

m 
.P_ 

2m 



(A2) 



A" 


a 




o 




o 


A 




— 





ig — T 

m 



(A3) 



For the respective gluon-antiquark vertices replace T — > T"' 
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Triple ultrasoft gluon - vertices^: 

AO 




, i,A 




10,C 




Heavy quark potential - vertices: 



(A4) 




Vr 



P 



P 



X 



p',fcf 



X 



-iV, 



(T) 



(P' - P) 



i Vc^'^ ® (p'-p)-[(^4-fc3)-(^2-fci) 



(p'-p)-* 



(fc4-fc3)^ + (fc2;~fcl)^ 
2{p'-pr 



2[(p'~p)(fc4-fc3)]^+2[(p'-p)(fc2^fcl)P 

(p'-p)e 



p, p'. fcf 

Heavy quark potential - ultrasoft gluon - vertex: 




[p' - P) • 
(P' - P)^ 



(A5) 



(A6) 
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